Self-trapped Exciton and Pranck-Condon Spectra Predicted in LaMnOa 
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Because the ground state has cooperative Jahn- Teller order, electronic excitations in LaMnOs 
are predicted to self-trap by local rearrangement of the lattice. The optical spectrum should show 
a Franck-Condon series, that is, a Gaussian envelope of vibrational sidebands. Existing data are 
reinterpreted in this way. The Raman spectrum is predicted to have strong multiphonon features. 
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In small molecules, excited electronic states generally have altered atomic coordinates, which leads to Franck- 
Condon multi-phonon sidebands in electronic spectra ^. In solids, electronic excited states are often delocalized, 
eliminating such effects H]. However, if excited states self-localize 0M| then Franck-Condon effects should reappear 
in the form of Gaussian broadening of the pure electronic transition B . Here we argue that such effects are crucial 
to a correct interpretation of LaMnOs, the parent compound of the "colossal magnetoresistance" materials. 

LaMnOa has a cubic to orthorhombic cooperative JT distortion M at 800K which corresponds to wavevector 
Q = (tt, tt, 0). This drives orbital order JTJ: x and y-oriented eg orbitals alternate in the x — y plane, in a layer structure 
shown schematically in Fig. |l|. A "minimal" model has two Mn d orbitals {eg states ipx^-y^ and ip^^^ -r'^) , Hubbard 
U , Hund energy J to align the Eg spin with the S — 3/2 spin of the 3 t2g electrons, a hopping term permitting band 
formation, and very important, electron-phonon coupling which allows local distortions of the oxygen environment to 
split the energies of the singly occupied eg levels and drive the JT transition. 

The highly-doped phases of LaMnOa show a remarkable interplay of charge, orbital, and spin order ||]. The relative 
importance of Coulomb, electron-phonon, and other effects is controversial. For light doping, however, charge order 
is not an issue and the problem simplifies. With no "empty" (Mn''+) sites for an eg electron to hop into, the large 
size of U suppresses hopping. The electron-phonon term becomes the dominant part of the remaining Hamiltonian, 
and is simple enough that low-temperature properties can be solved. We have already reported |S[| properties of the 
self-trapped hole ( "anti Jahn- Teller" small polaron) which forms at low doping. Its sensitivity to spin order influences 
the magnetic states seen at low doping. We also gave a preliminary treatment of electronic excitations |1^. Here we 
give a more complete solution for the lowest electronic excitation, finding a self-trapped exciton, insensitive to spin 
order, whose Franck-Condon effect should dominate cr(i^) and the resonant Raman spectrum, and might also appear 
in luminescence. These results follow cleanly from our model, but have not previously been recognized as relevant. 
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FIG. 1. a. schematic of Jahn- Teller distorted LaMnOa x — y plane with orbitals ordered, b. The lowest electronic excitation 
with an orbital rotated and the lattice relaxed. 

We make the strong-coupling approximation U —* oo hy considering only states with no double occupancy of eg 
orbitals; this suppresses hopping and leaves only two terms in the Hamiltonian. Vibrational degrees of freedom are 
modelled by allowing each oxygen to move along the nearest-neighbor Mn-O-Mn axis. 
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Here ui^x is the displacement from cubic perovskite position of the nearest oxygen in the x-direction to the Mn 
atom at L The Jahn- Teller energy is modelled by a linear reduction of energy of an occupied 3x^ — r^ orbital if the 



corresponding two oxygens in the ztx direction expand outwards, and similarly for y and z oxygens if 3y^ — r^ or 
3z^ — r^ orbitals are occupied, 

Ti-JT = -g'^fli^aiui^a - Ui^^a)- (2) 

Here fii^x is the occupation number of the 3x^ — r^ orbital, and U£^^x is the displacement of the oxygen nearest the 
Mn atom at i in the — x-direction. Since no hopping occurs in this model, the type and degree of magnetic order is 
irrelevant, and no change occurs in the electronic spectrum a s T is increased through the Neel temperature Tjv=140K. 
The adiabatic parameter is a = Hwq/A where loq = \J KjM K,lh meV is the oxygen vibration frequency lO] 
and 2A «1.9 eV is the Jahn- Teller gap (explained below.) Since a is small, we solve the Hamiltonian in adiabatic 
approximation, then add quantized lattice vibrations. There is an infinite set of equally good distorted ground states 
[ jl2[ whose degeneracy is broken by anharmonic terms left out of the minimal H. — Tivih + 'Hjt . We simply adopt the 
distortion seen experimentally and shown in Fig. |l|. The A sublattice is defined to be Mn sites where exp{iQ • £) is 1, 
and the B sublattice to be sites where it is -1. The resulting ground state electronic wavefunction M is 

A B 

\o>=Y[cUe)Y[^Unm>, (3) 

where |{0} > refers to the lattice state with oxygens displaced by ±uo = '^g/K as shown in Fig. [^ and in their 
vibrational ground states. The orbitals 4'x and -0^ which are occupied on the A and B sublattices are orthogonal x- 
and y-directed states ipxy = {—ipsz^-r^ i ip^2_y2)/\/2. In state (|^) each occupied Mn orbital has its energy lowered 
by A = 8g^/K and each unoccupied orbital has its energy raised by A. This costs elastic energy Kuq = A/2 per Mn 
atom (two of the three oxygen atoms per Mn site are displaced) giving a total condensation energy of —A/2 per Mn. 



Our picture is substantially the same as the earlier theory by Millis |12|. 

If oxygen atoms were frozen in these optimal positions, the lowest electronic excitation (an "orbital defect," Fig. Ill 
b) would cost 2 A = 1.9 eV, the JT gap. This would form a narrow "orbiton" band ||l^ because terms of order t^ /U 
(neglected here) allow the orbital defect to exchange sites, lowering the energy by r~/U. However, there is a much 
greater energy lowering if instead the lattice locally undistorts, pinning the orbital excitation on a single site as 
shown in Fig. Qd. This reduces the energy of the orbital defect from 2A to A. Thus the lowest electronic excitation 
is a strongly self-trapped exciton, somewhat like the Frenkel exciton seen in molecular crystals [0. 

If this excitation is excited optically, the Franck-Condon principle applies, and a sequence of vibrational sidebands 
of the localized exciton will appear in the spectrum. The optical conductivity is 

cr(w) = — t: — 22 7 1 < •^"'1^ ■ ^*"' > '^"^(^ + ("-/ ~ n,)huJo - ?iw), (4) 

n,n' 

where il is the volume of the crystal, and the number of cells N appears because the localized excitation can be 
created on any Mn atom. The partition function Z is {1 — exp{f3huJo))~^ because 6 oscillators couple to each electronic 
transition. The required dipole matrix element is 

<fn'\i-r^in>^ j (fR j dfxn'{R- RfW{r,R)eL ■ ri^x{r,R)Xn{R~ Ri)- (5) 

The crystal starts in state \in > with an electron at site in the electronic ground state \i >~ \X > with equilibrium 
coordinates Ri and vibrational quanta of the 6 surrounding oxygens denoted by {n} = (rii, n2, ns, 714, 715, ng). It ends 
up at the same site in final state \fn' > with the electron in state |/ >= \Y > of energy A with new equilibrium 
positions Rf and new vibrational quanta {n'} = (n'j^, 71,2, 71.3, 71.4, 71,5, rig). The notation rii denotes the total initial 
vibrational level tii + . . . + ng, and Uf is the total final vibrational level. The electronic part of the matrix element 
is a d to d transition which is forbidden when the surroundings are symmetric (as in the starting coordinates Ri or 
the final coordinates Rf if 77^ — nf = 0) . The origin of the observed optical transitions has been discussed by various 
authors p3Hl5[. In our localized picture with lattice displacements as primary factors, it is natural to recognize that 
any asymmetric oxygen breathing displacement will cause the Mn eg orbitals to acquire an admixture of 4p character. 
A typical mixing coefficient is 

f dV 

7== / dfV'3^2_^2— V'z/(e<i-ep) (6) 



where ipz is an orbital of p-character, partly Mn 4p, and partly oxygen 2p; dV/du^ is the perturbation caused by a 
displacement of oxygen 5 in the z-direction. The corresponding allowed optical matrix element is 



dr'4)^z2_^.2zi}z. 



(7) 



The resulting matrix element, Eqn. [g, is 

"yd 
< fn'\e ■ r]in >= — < n[n'2n'^n'^n'^nQ\[{ui + U2)ex + (us + U4)ey - 2{u^ + M6)e^]|"-l"2'^■3«4?^5"6 > (8) 

The general term (^ is complicated; we have evaluated only when the initial state has rii = or 1 quanta of 
vibration. Two or more quanta occur at <6% probability at 300K. Our approximate optical conductivity is then 
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The vibrational overlap integrals needed for cro are 
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and similar for states 2,3,4 except that for states 2 and 4, ~uq is replaced by +uq; this changes the sign of the prefactor 
of Eqn. O. We find that for both ctq and cri, (t^^ — ayy = a^zj A. The answer for ctq is 



" n=0 



-^)^(t + ("+^)"° 



(12) 



The formula for ui is more complicated and will not be given here. 

Our theory compares reasonably well with experiment. (j{ijj) has been measured by reflectivity on polycrystalline 
samples at room temperature |0,nq|, single crystals at low T p7| ], and cleaved single crystals at room temperature 
|18| , all with consistent results. Jung et al. |16| have expressed their cr(w) as a sum of Lorentzian peaks. The lowest 
peak is centered around 1.9 eV with a width of 1 eV. Fig. shows their Lorentzian fit, and compares with our theory, 
with 2A = 1.9 eV and amplitude chosen to agree with experiment. 
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FIG. 2. Optical conductivity of LaMnOs. The points are the lowest Lorentzian oscillator fit by Jung et al. to their data. 
The dashed curve is a T = sum of convolved Lorentzians centered at the vibrational replicas shown as vertical bars; the solid 
curves are T = (lower) and T — 300K (upper) sums of convolved Gaussians, also shown in the inset on a logarithmic scale. 
Tick marks in the inset denote decades. 



The intensity seen experimentally |lj] corresponds to /=0.16 oscillators per Mn atom in the lowest peak, where 
J dwa{oj) is defined as {irNe^ /2mil)f, with m the electron mass. Our theory gives / = 2{'-fd)'^{m/M)(2A/hujQ + 1) 
and has the value 0.18{'~fd)'^ x 10^^. To agree with experiment, jd must be 10. This is reasonable since the electron- 
phonon mixing coefficient may be a few per A oxygen displacement, while the dipole matrix element is likely to be 



a few A. Most authors p3-17| agree in assigning the 2 eV structure to eg to Cg transitions activated by mixing of 
p-character, usually from surrounding oxygens. However, the origin of the mixing is usually taken as the dispersive 
broadening of the Cg bands; in our strong-coupling picture, these effects enter at order t/U and are neglected, but 
might be comparable to the phonon effects which we calculate explicitly. 

The delta function peaks of our model theory should be replaced by local densities of vibrational states on oxygen 
atoms. If this is mimicked by using a sequence of convolved Lorentzians {L, L * L, . . .) the result agrees closely 
with experiment. A sequence of convolved Gaussians may be more realistic, and is shown in Fig. as a dashed 
line and enlarged on a logarithmic scale. Our theory is not the only candidate; band theory |1^] yields reasonable 
agreement at T = 0. However, our theory makes definite predictions which differ from those of band theory. We 
find the onset of absorption at T = to be the one-phonon structure at A + cjph ~1 eV, weaker by 10'* than 
the 2A=1.9 eV peak absorption. The fine structure of subsequent phonon peaks will be harder to resolve because 
of increasing multiplicity of vibrational quanta available for multi-phonon absorption. Another definite prediction 
concerns temperature dependence. Contrary to approaches using dispersive bands [[l9| , |l5[ , <t{llj) in our approach is 
not affected by the loss of magnetic order as temperature increases. Contrary to band theory, we predict very weak 
temperature effects near the absorption peak at 2A, consistent with experiment [ |l7|Jl8| . By contrast we predict big 
changes of intensity in the weak features near w = A, including new weak peaks at A — nwph for n = 0, 1, . . . activated 
by temperature, and shown in the inset to Fig. g. These effects are too small to be seen in existing reflectivity data, 
but should be measurable by absorption in very clean thin films. If electronic exitations survive long enough to relax 
by luminscence, then a very large Stokes shift is predicted. 

Another manifestation of the Franck-Condon effect should appear in the Raman spectrum pffl, particularly if the 
laser frequency is near the 2A=1.9 eV Jahn- Teller gap edge. Multi-phonon Raman features should appear with 
intensity similar to the one-phonon spectrum. The incident photon creates (among various other virtual excitations) 
a self-trapped exciton in a superposition of multi-phonon states (as in the left inset of Fig. 0.) The virtual exciton 
can reemit a photon, returning either to the ground state or (with nearly equal amplitude) to various one-phonon 
or multi-phonon excited states. The theory was given by Shorygin pl| . A particularly clear solid-state example is 
Martin's |g^ observation of multiples of the localized phonons of the Mn04^ impurity complex in Csl. Multiples 
of the Jahn- Teller related vibrations should appear in the Raman spectrum of LaMnOa, but less distinct than in 
Martin's work since the local phonon density of Jahn- Teller vibrations will be less localized and correspondingly 
less well-structured. Published Raman measurements ||ll| on undoped LaMnOs do not extend to the multi-phonon 
region. An unpublished two-phonon replica of the Jahn- Teller phonon in LaSr2Mn207 is seen by Romero et al. [ p3| , 
and similar features are seen in unpublished p4| work on LaMnOa. These are probably the effect we are predicting. 
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